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ABSTRACT 


Based on new ISO mid-infrared observations and ground based Ha and 
near-infrared observations, we report the detection of a bright starburst in the 
intragroup medium (IGM) of the famous compact group of galaxies Stephan’s 
Quintet (Source A in Fig.l). We demonstrate that this starburst is caused 
by a collision between a high velocity ((5V~ 1000 km/sec) intruder galaxy 
(NGC7318b) and the IGM of the group. While this is the only starburst 
known today that is induced by a galaxy/cold-intergalactic-medium collision, it 
provides new constraints to the theory for interaction-induced starbursts, and 
may hint at a new mechanism for the star formation excess seen in more distant 
clusters. 

Subject headings: galaxies: interactions - galaxies: intergalactic medium - 
galaxies: ISM - galaxies: starburst ~ galaxies: active - infrared: galaxies - stars: 
formation 



1. Introduction 


Starbursts (Rieke et al. 1980; Gehrz et al 1983) are “cosmic fireworks” where hundreds 
of thousands of massive stars that are up to a million times brighter than the Sun are born 
and die in a short time scale of a few 10^ years (compared to the age of most galaxies of 
~ 10^'^ years). Despite the wide interest, many basic questions about starbursts remain 
unanswered. First of all, how are they triggered? Theoretically the only well established 
triggering mechanism for starbursts involves low velocity enconnters/mergers between 
galaxies (Larson & Tinsley 1978; Joseph et al. 1984; Lonsdale et al. 1984; Noguchi & 
Ishibashi 1986; Mihos & Hernquist 1996). In this letter, we shall show that a bright 
starburst in the intragronp medinm (IGM) of the famous compact group of galaxies 
Stephan’s Qnintet is triggered by a new mechanism, namely a collision between a high 
velocity ((5V~ 1000 km/sec) intruder galaxy (NGG7318b) and the IGM of the gronp. 

Stephan’s Quintet (hereafter SQ), discovered 120 years ago (Stephan 1877), is perhaps 
the most famous example of the compact galaxy group class that is characterized by 
aggregates of 4-8 galaxies with implied space densities as high as those in cluster cores 
(Hickson 1982). SQ has been observed with almost every large telescope and in almost 
every domain of the electromagnetic spectrum. These data suggest (Moles et al. 1997) that 
SQ has been visited by at least two “intrnders” (NGG7320c and NGG7318b) in the past few 
10® years. Gonsequently, its intragronp medinm (IGM) is hlled with HI gas stripped from 
its spiral members (Shostak et al. 1984). In an attempt to study the star formation activity 
and the dust properties of this IGM we made observations at 11.4/im, 15/im (ISOGAM), 
60/im and lOO/rm (ISOPHT) using the Infrared Space Observatory (Kessler et al. 1996). In 
these observations, an outstanding bright source (Source A in Fig.l) is detected in a region 
qnite far away from centers of member galaxies. In this letter, we pnblish the ISOGAM 
15/im data (the rest of ISO data will be published in subsequent papers) which, together 



with data from the following-up gound based Hq, and NIR observations, demonstrate 
conclusively that this bright source is associated with a starburst in the IGM of SQ. 


2. Observations 

2.1. ISOCAM observations 

The ISOCAM observations at IS/rm (Fig.l) were made on 1996 May 23 using the 
ISOCAM Long-Wavelength-Channel array (32x32 pixels). Raster scans with PFOV = 

M = 12, 5M = 48'' (in-scan) and N = 20, 5N = 6" (cross-scan) were made in the 
micro-scan mode (CAMOl) using hlter LW3 (Aq = IS/rm, 5\ = 6/im). The la rms noise 
measured in the inner part (~ 5' x 3') of the image is 0.04 mJy/pixel, and is higher near the 
edges of the image. The absolute calibration is taken from the ISOCAM Observer’s Manual 
(1994), which has an uncertainty of less than 30% (Sauvage, private communication). The 
angular resolution (FWHM) of the image is 10". The basic data reduction was done using 
the CAM Interactive Analysis (CIA) software^. 

2.2. Ha observations 

The Ha observations (Fig.2) were made in August 1997 using the 2.2m telescope of 
the observatory of Max-Planck-Inst it ut at Calar Alto (Spain). In order to separate the 
Ha emission associated with the intruder galaxy NCC7318b (n ~ 5700 km/sec) from the 
Ha emission associated with the rest of the group (n ~ 6600 km/sec), CCD images of the 
Ha-[NH] emission were obtained using two narrow band hlters 667/7 and 674/8, centered at 

^CIA is a joint development by ESA Astrophysics Division and the ISOCAM Consortium 
led by the ISOCAM PI, C. Cesarsky, Direction des Sciences de la Matiere, C.E.A., France. 
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6667A and 6737A with FWHM bandwidths of 66A and 76A, respectively. The transmission 
of the 667/7 hlter for the Hq, line of redshift 0.019 {v = 5700 km/sec) is 0.49, and for the 
Hq, line of redshift 0.022 (6600 km/sec) is 0.04. The transmissions of the 674/8 hlter for the 
same two components of the emission are 0.06 and 0.49, respectively. The 5700 km/sec 
component of the Hq,-[NII] emission is therefore estimated from the image obtained using 
the 667/7 hlter, and the 6600 km/sec component from the 674/8 image. The continuum 
in both images is subtracted using an R-band CCD image of the same held. To enhance 
the sensitivity for the dihuse emission, both Hq,-[NII] maps (original resolution ~ 1") are 
smoothed to a round Gaussian beam with FWHM=2". Note that in some regions the two 
components have similar morphology. This could be an indication that the emission in those 
regions is rather associated with another component, namely the one with recession velocity 
of 6000 km/sec (Moles et ah 1997; Shostack et al. 1984). However this uncertainty does 
not ahect seriously the Hq, huxes (which are sums of the 5700 km/sec and the 6600km/sec 
components) reported in Table 1, because while the 6000km/sec component might have 
been counted twice, its flux would have been underestimated by about a factor of two in 
both the 667/7 map and the 674/8 map, given that the transmissions of the two hlters for 
the 6000 km/sec component are only about half of the transmissions for the 5700 km/sec 
and the 6600 km/sec component, respectively. 


2.3. Near Infrared K’-band observation 

A near infrared K’-band (Aq = 2.1pm) image of the central part of Stephan’s Quintet 
were obtained using the 3.5m telescope of the observatory of Max-Planck-Institut at Calar 
Alto (Spain), mounted with the 256x256 pixels (pixel size =0."81) MAGIC camera. The 
observations were made in March 1997. The weather conditions were photometric and 
stellar photometry consistent to ~ 2% throughout the night. The seeing was about 1". The 
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image covers about 4' x 4' sky (Fig.2). In the central 3' x 3' field, the la rms noise is 20.4 
mag/arcsec^. Because of the mosaic, the noise is higher near the edges. 


3. Results 

3.1. A starburst (Source A) far away from galaxy centers 

Figure 1 presents the central ~ 4' x 4' region of the 15/im ISOCAM image of 
Stephan’s Quintet (false color). As shown by Sauvage et al. (1996), the IS/xm emission is 
a good, nearly extinction free star formation rate indicator. Located slightly upper-right 
from the map center. Source A lies about 1' (~ 25 kpc at H = 80 Mpc where Ho=75) 
north of NGC7318b and about 1' west of NGC7319. It covers a region of ~ 30" and its 
peak surface brightness in the 15/im map is second only to that of the Sy2 nucleus of 
NGC7319. At 80Mpc the 15/im flux density of Source A corresponds to a luminosity 
Li 5 ^ = 4.8 IO^Lq, about 30 times higher than the (1.7 IO’^Lq) of the foreground Sd 
galaxy NGC7320 (assumed to be at lOMpc; Table 1). 

In an early photographic observation including both the Ha emission and the underlying 
continuum emission, Arp (1973) identified several HII regions in the region of Source A. Hq, 
emission was also detected in the Source A region in a recent published Hq, image centered 
on the Sy2 galaxy NGC7319 (Moles et al. 1997) though Source A falls on the edge of that 
image and therefore is not fnlly covered. A grid photometry map for SQ by Schombert et 
al. (1990) reveals a concentration of very bine data points (B—V= 0.3-0.5) in this region. 
The Hq knots and blue color provide strong support for an interpretation that the MIR 
emission is due to dust heated by massive young stars. 

In order to firmly establish the starburst nature of Source A, we carried out new Hq 


and K' (2.1/im) imaging observations using an optical GGD and a near-infrared detector 
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array, respectively (Section 2.2 and 2.3). With two narrow band filters (667/7 and 674/8), 
the 5700 km/sec component and the 6600 km/sec component of the Hq, emission, the former 
being associated with the current intruder galaxy NGC7318b and the latter with the rest 
of the group (Table 1), are separated. In Fig.2 the contours of the two Hq, components are 
overlaid on the K' image. Indeed strong Hq, emission from both components is detected 
in the Source A region, and the position of the Hq peak (associated with the 6600km/sec 
component) coincides well with that of the 15/rm emission peak. The Hq equivalent line 
width, EW(Hq), of the region is 115±15A. Such high EW(Hq) is found only in starbursts 
(Keel et ah 1985). 

The star formation rate (SFR) in the Source A region can be estimated from both 
the 15/im and the Hq luminosities. A ratio between flux densities of the 15/im and of the 
thermal radio at 6cm fi 5 ^/fth, 6 cm = 85 (mJy/mJy) is assumed, taken from the median value 
of the star formation regions in M51 (Fig.2b of Sauvage et al. (1996)). Following Kennicutt 
et al. (1994) and adopting a Salpeter IMF with mi = 0.1 M© and = 100 M©, we obtain 
from Source A’s 15/rm luminosity (4.8 10® Lq) an SFR of 0.81 MQ/yr. Assuming the same 
IMF, the Hq luminosity (2.2 10^ Lq) gives an SFR of 0.66 MQ/yr. The two estimates of the 
SFR agree with each other within the uncertainties. This also suggests that the extinction 
of the Hq emission in Source A is low. 

For starbursts, according to Bruzual & Chariot (1993), the stellar-mass to NIR 
luminosity ratio is insensitive to the starburst age, and therefore the NIR luminosity is a 
good stellar mass indicator. Using the model of Bruzual & Chariot (1993) and assuming 
that the stellar population in Source A is the product of a starburst of age 10^ — 10® years 
(the widely accepted age span of starbursts), the stellar mass in Source A can be estimated 
from its K' magnitude (Table 1), which is in the range 0.8 — 1.6 10^ Mq. Dividing this 
mass by the star formation rate derived above, one hnds that the starburst in Source A is 



only ~ 1—2 10^ years old. Any possible contribution from an underlying older population 
(e.g. from stripped stars) may only reduce the estimate of this age, because less NIR light 
would then be due to stars formed in the current burst. 


3.2. Triggering mechanism: IGM—intrnder collision 

How is the starburst in Source A, which is more than 20 kpc away from centers of 
nearby galaxies, triggered? HI observations of SQ (Shostak et al. 1984) reveal three velocity 
components associated with the group (not including the foreground galaxy NGC7320): 1) 
a 6600km/sec component distributed throughout the intragroup space, 2) a 6000km/sec 
component north of N7318b and 3) a 5700km/sec component located south of N7318b. 
Moles et al. (1997) suggested that both the two latter components are associated with 
NGC7318b, the velocity difference between them due to the rotation. Both the 6600 and 
6000 km/sec components show local maxima in the region of Source A, both with HI 
surface density of ~ 3.3 10^° atoms/cm^. The X-ray (Pietsch et al. 1997) and the radio 
continuum (van der Hulst & Rots 1981) observations show a shock front, generated by the 
collision between NGG7318b and the IGM, on the eastern side of NGG7318b (also visible 
on our Ha maps, Fig.2). Noticeably, Source A is located in a region without detected X-ray 
emission. ROSAT maps (Pietsch et al. 1997) show that Source A lies either: 1) in a gap in 
the shock related X-ray emission (if the detached emission NW of Source A is interpreted 
as an extension of the shock front) or 2) on the N edge of the shock front. 

The HI, X-ray and radio continuum data in the literature and our MIR, H^ and 
NIR data are all directing to a scenario such as the following: During the intrusion of 
NGG7318b into SQ, a high speed (~ 1000 km/sec) collision between the local maxima of 
the two HI components in Source A region, one (6600 km/sec component) associated with 
the intragroup gas and the other (6000 km/sec component) with the intruder NGG7318b, 
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occurred ~ 10^ years ago. This collision is perhaps part of the larger scale collision which 
has resulted in the shock front seen in the X-ray and radio continuum observations on the 
eastern side of NGC7318b. The collision triggers an instantaneous starburst, giving rise to 
the strong MIR and Hq, emission (its is about twice of that of the starburst in the 
collision region of the Antennae Galaxies (Vigroux et ah 1996)). The hot gas produced by 
the collision has already cooled in the region of Source A presumably due in part to the 
locally higher gas density which results in a higher cooling rate. 

Interestingly, Source A sits also in a region where two long and faint optical arms 
seemingly run across with each other (Arp & Lorre 1976). Both arms are likely being 
stretched out from NGG7318b, with the NGG7318a (an early type galaxy, Table 1) in 
the background/foreground of the western arm (Moles et al. 1997). It is tempting to 
assume that these arms are products of tidal interactions between NGG7318b and the 
other two galaxies (NGG7318a and NGG7319), and the starburst in Source A is triggered 
by the interaction of these two arms of NGG7318b. Strong evidence against this scenario, 
and favoring the IGM-intruder collision scenario, is the presence of a strong 6600km/sec 
component, which is in fact the dominant component with apparent connection with the 
shock front, in the Hq emission in Source A (Fig.2). If there is no participation of the IGM 
in the starburst, the 6600 km/sec component is not at all expected. The fact that the Hq 
emission in the Source A region is dominated by the 6600km/sec component conclusively 
demonstrate that most of the star formation in Source A region must be happening in the 
IGM. 

In addition, the time scale for the development of tidal tails is more than a few 10® years 
(Barnes & Hernquist 1992; Elmegreen et al. 1993) while the crossing time of NGG7318b is 
only a few 10'^ years, too short for the tidal effects to act. Indeed it is more likely that the 
two optical ‘arms’ are not ‘tidal arms’ but are, together with Source A, products of some 
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hydrodynamical processes induced by the high speed collision between the intruder and the 
rest of the group. The eastern ‘arm’ is almost certain to be associated with the shock front 
(Moles et ah 1997) although the mechanism for the western ‘arm’ is still not very clear. 

Once the IGM environment of the starburst in Source A is determined, the next 
question is: could it be due to star formation in the cooled gas in the aftermath of the high 
speed collision, similar to that proposed for the star formation in ‘cooling flows’ (Fabian 
et al. 1991)? It has been strongly argued, because of the red colors of the ‘cooling flow’ 
galaxies, that if star formation does occur in ‘cooling flows’ it has to be severely deficient in 
massive stars (Fabian et al. 1991). This is contrary to what we are seeing in Source A where 
the high EW(Hq,) clearly indicates a starburst with a lot of massive stars. Hence the ‘star 
formation in cooling flow’ mechanism is irrelevant here. It is more likely that the starburst 
in Source A is triggered by the gravitational instability casued by some hydrodynamical 
processes which in turn are induced by the high velocity collision between the IGM and the 
intruder. This might have some similatities, albeit in a much larger scale, to the triggering 
mechanism of the so called ‘self-propagating’ star formation mode (Blaauw 1964) in the 
Galactic disk, which is associated with the shocks generated by HII regions (Elmegreen & 
Lada 1977) and by supernova remnants (Herbst & Assousa 1977). 

In summary, our conclusion that the starburst in the Source A region is triggered by 
the high speed (~ 1000 km/sec) collision between the IGM and the intruder NGG7318b is 
supported by the following facts: 

1. The Hq, data show that the star formation in Source A is occurring both in the IGM 
(the 6600 km/sec component) and in the intruder (the 5700 km/sec component). In 
particular, the fact that the Hq, emission in the Source A region is dominated by 
the 6600km/sec component conclusively rules out any interpretation involving only 
processes conhned within the intruder NGG7318b. 
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2. The age of the starburst 1— 2 10'^ yrs) estimated from the MIR, and NIR data 

is consistent with the dynamical time scale of the high speed collision ( 10^ yrs). 

The probability that both the burst and the collision are happening simultaneously 
(within such a short time scale) would be very low if the former is not causally related 
to the latter. 

3. The HI observations of Shostak et ah (1984) show that both the intragroup gas and 
the gas associated with NGC7318b have local concentrations in the Source A region. 
These cold gas components are not only the participants of the high speed collision 
which triggers the starburst, but also the reservoirs of the necessary material for the 
starburst to proceed. 

4. The morphology of the 6600km/sec component of the emission shows a possible 
connection between Source A and the shock front, indicating that Source A is likely 
part of a larger scale collision that is still going on between the IGM and the intruder. 

5. Although Source A has no detectable X-ray emission, it is surrounded by X-ray 
emitting hot gas, again indicating a connection between the shock front (with which 
the X-ray emitting hot gas is associated) and the starburst. The high pressure 
imposed by the surrounding hot gas may have also helped sustain the starburst by 
conhning the cold gas in that region. 


4. Discussion 

4.1. Theoretical implication 

Theoretical works on interaction-induced starbursts are still in their early stages 
(Mihos & Hernquist 1996). Most of the studies have concentrated on the effects of the 
gravitational tidal force, while the effects of the hydrodynamic shocks generated by the 
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collision of two gas condensations have been left out (Noguchi & Ishibashi 1986; Mihos et 
ah 1993). The scenario envisionized by Jog & Solomon (1992) for starbursts induced by 
low velocity (~ 100—300 km/sec) galaxy-galaxy collisions is difficult to apply to Source A. 
In that scenario a starburst occurs when the Giant Molecular Clouds (GMC’s) preexisting 
in the two colliding disks are compressed by the hot shocked gas. In the case of Source 
A, it is not clear whether the IGM contains GMC’s (the CO survey of SQ by Yun et 
ah (1997) unfortunately did not cover the Source A region). Even if there are indeed 
preexisting GMC’s in the Source A region, the colliding velocity is so high (~ 600—1000 
km/sec) and close to the typical escape velocity in galaxies that the shocked gas may be 
left behind when the colliding parties move away from each other while GMC’s are not 
affected by the collision because of the low hlling factor (Jog & Solomon 1992). Clearly 
more theoretical work, in particular simulations with a full treatment of gas dynamics, is 
needed to understand what is going on in starbursts such as Source A. 


4.2. A new mechanism for Butcher—Oemler effect? 

Events like Source A are very rare in the local universe. However in the earlier universe, 
in particular in young clusters where high velocity encounters involving late type galaxies 
are frequent and a large part of the intracluster space may be hlled by cold stripped gas, 
collision between gas rich galaxies and the cold intracluster medium (ICM) may be quite 
common. Whether the starbursts so induced are at least partially responsible for the 
enhanced star formation in those young, distant clusters (Butcher & Oemler 1984), needs 
to be further explored (for an alternative mechanism, see Moore et al. 1996). 
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Table 1. 



R.A. 

Dec. 

Type 

Vr 

f 1 


EW(H«)' 

(J2000) 

(J2000) 


km/sec 

mJy 

10 -13 g 2 ;g/cm^ /sec 

A 

NGC7319 

22h36m03.5s 

+33d58'33" 

SBb (Sy2) 

6764 

79.8 

1.41T0.30 

7.6±2.5 

NGG7318a 

22h35m56.7s 

+33d57'56" 

Epee 

6630 




NGG7318b 

22h35m58.4s 

+33d5r58" 

SBbe 

5774 




NGG7318a/b5 





19.3 

1.08±0.16 

5.0T1.2 

NGG7317 

22h35m52.0s 

+33d56'41" 

E 

6599 

2.1 



NGG7320 

22h36m03.5s 

+33d56'54" 

Sd 

786 

27.5 



Source A 

22h35m58.7s 

+33d58'55" 

SF region 


11.9 

1.27±0.13 

115±15 

Source B® 

22h36ml0.2s 

+33d57'21" 

SF region 


1.8 

0.14±0.01 

85±10 


^The 15/im flux density, measured through aperture photometry covering areas included in 2a contor 
uncertainties are about 30%, dominated by the errors in the calibration. 

flux, corrected for the [NIIA6583/A6548] contamination. For NGC7319, the Sy2 galaxy, t’ 
[NIIA6583 ]/Hq = 1.8 is taken from Keel et ah (1985). For other galaxies and star formation reg 
[NIIA6583 ]/Hq ratio is assumed to be 0.4 (Kennicutt et ah 1994). The errors are dominated by the unce 
in the continuum subtraction. No extinction correction is applied. 
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Fig. 1.— ISOCAM LW3 (15/im) image (false color) of Stephan’s Quintet field. The scale is in 
logarithm and the units are 0.1 mJy/pixel (pixel= 6" x 6"). As shown by the chromatometer, 
the range of the surface brightness shown in the hgure is —1 < log[S'i 5 ^m/(0.1 mJy/pixel)] < 
2. The hve members of the group, including the foreground galaxy NGC7320, are clearly 
detected. Four other sources outside the galaxies are marked in the image: Source A and 
Source B are star formation regions associated with the group (see Table 1). Source C is a 
foreground star, and Source D is a background galaxy. 


Fig. 2.— K' (2.1/im) image overlaid by contours of the Hq,-[NII] emission, covering the same 
held of sky as in Fig.l. The K' image is in logarithmic scale. The two components of the Hq,- 
[NII] emission are marked by different colors: the contours of the 6600 km/sec component are 
in red and those of the 5700 km/sec component in blue. The two contour sets have the same 
levels, which are 0.3, 0.9, 2.7, 8.1, 24.3 10“^® erg/cm^/sec/pixel (pixel= 0.533" x 0.533"). 
The hrst contour level is chosen to hide most of the residuals of the continuum subtraction, 
the most important source of the errors of the Hq,-[NII] maps. 
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